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This paper describes the results of an analysis of localized catalytic heating effects to 
the U.S. Space Shuttle Orbiter Thermal Protection System (TPS). The analysis applies 
to the High-temperature Reusable Surface Insulation (HRSI) on the lower fuselage and 
wing acreage, as well as the critical Reinforced Carbon-Carbon on the nose cap, chin panel 
and the wing leading edge. The object of the analysis was to use a modified two-layer ap- 
proach to predict the catalytic heating effects on the Orbiter windward HRSI tile acreage, 
nose cap, and wing leading edge assuming localized highly catalytic or fully catalytic sur- 
faces. The method incorporated the Boundary Layer Integral Matrix Procedure— Kinetic 
(BLIMPK) code with streamline inputs from viscous Navier-Stokes solutions to produce 
heating rates for localized fully catalytic and highly catalytic surfaces as well as for nomi- 
nal partially catalytic surfaces (either Reinforced Carbon-Carbon or Reaction Cured Glass) 
with temperature-dependent recombination coefficients. The highly catalytic heating re- 
sults showed very good correlation with Orbiter Experiments STS-2, -3, and -5 centerline 
and STS-5 wing flight data for the HRSI tiles. Recommended catalytic heating factors were 
generated for use in future Shuttle missions in the event of quick-time analysis of damaged 
or repaired TPS areas during atmospheric reentry. The catalytic factors are presented 
along the streamlines as well as a function of stagnation enthalpy so they can be used for 
arbitrary trajectories. 
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Dij 

multicomponent diffusion coefficient 

h° 

enthalpy of formation 

H w 

wall enthalpy 

hare 

heating parameter 

H C 1 

arc-jet nozzle centerline enthalpy 

Mi 

normal mass flux of i th species 

p 

- 1 meas 

measured test article surface pressure 

q 

heat flux per unit area 

T 

temperature 

t 

time 

(Xo,Yo,Zo) 

axial, lateral, and vertical Shuttle coordinates 

N 

atomic nitrogen 

n 2 

molecular nitrogen 

NO 

nitric oxide 

0 

atomic oxygen 

0 2 

molecular oxygen 

x/L 

Orbiter axial distance from apex normalized by reference length (L=1290 inches) 


I. Introduction 

T he capability to analyze the catalytic heating effects of localized repaired or damaged surfaces on the 
Thermal Protection System (TPS) of the U.S. Space Shuttle Orbiter was a major recommendation 
set forth by the Columbia Accident Investigation Board 1 * before the Orbiter could be returned to flight. 
Because exposed underlying layers of material resulting from damage to the TPS surface can have different 
thermochemical properties from the Reaction Cured Glass (RCG) coating on the Orbiter windward acreage 
and/or the Reinforced Carbon-Carbon (RCC) on the wing leading edge, the chemical surface heating, or 
catalytic heating, occurring during Orbiter atmospheric reentry can significantly increase in the localized 
damage area beyond the values predicted for flow in chemical equilibrium. In addition, any repair to the TPS 
may cause a change in the chemical make-up of the repair site, which can also cause higher catalytic heating 
compared to the surrounding TPS. In order to incorporate localized catalytic heating into the Shuttle TPS 
thermal math models, catalytic heating amplification (or bump) factors will be used to increase the nominal 
heating rates to damaged or repaired areas. 

In support of Return to Flight (RTF), the National Aeronautics and Space Agency (NASA) at Johnson 
Space Center (JSC) in Houston, TX requested that analyses be performed to assess the localized heating 
effects from highly catalytic surfaces that may be present on the Orbiter TPS during atmospheric re-entry. In 
the event that the Orbiter becomes damaged and must re-enter the atmosphere with TPS damage exposure 
or repair sites, a quick surface heating analysis must be performed to determine if aerothermal heating will 
become problematic. Because of the lack of catalytic material test data and because performing parametric 
analyses of catalytic heating solutions in localized areas with computational fluid dynamics (CFD) models 
is time and resource intensive, a modified two-layer approach consisting of streamline inputs derived from 
chemical non-equilibrium Navier-Stokes solutions and a boundary layer method was used to perform the 
analyses. The analyses were compared to the limited amount of catalytic flight data obtained from the 
Orbiter Experiments (OEX) program. 2, 3,4 

The objective of the analyses was to ascertain whether the thermo-chemical kinetic gas model contained 
in the Boundary Layer Integral Matrix Procedure - Kinetic (BLIMPK) code 5 could be used to predict the 
catalytic heating rate overshoots exhibited in the OEX program flight measurements with temperature- 
dependent recombination coefficients from material data analyses of RCG 6 and C742 7 coated tile. Early 
BLIMPK analyses used parametrically adjusted constant recombination rates to match the radiometer flight 
data at the Orbiter stagnation point. 8 Later analyses for the Shuttle Entry Air Data System (SEADS) nose 
cap revealed that localized catalytic surfaces could produce heating rates that overshoot the theoretical fully 
catalytic (equilibrium) values. 9 Additional observations of this overshoot phenomenon were also documented 
in other Shuttle analyses. 10 If the modified two-layer method could be validated, it would be used to assess 
surface catalytic heating effects to damaged or repaired TPS locations on the Orbiter during atmospheric 
re-entry in support of RTF recommendations assuming a localized fully catalytic damage site. 
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II. Methodology 


The surface catalytic heating study presented in this paper was conducted over approximately two years 
and included research of historical flight data and analyses, arc-jet tests on coated and uncoated High- 
temperature Reusable Surface Insulation (HRSI) tiles and RCC materials, a study of analytical approaches 
and a large number of new analyses with modified analytical methods and procedures. The following sections 
describe the progression of the catalytic heating study and the methods employed to develop a working 
analytical model for obtaining localized catalytic heating factors on a damaged or repaired Orbiter heat 
shield. 

A. Analysis of Flight Data 

The first study of catalytic heating factors for 
use on damaged or repaired locations on the Shut- 
tle TPS began in late 2003. Because no measure- 
ments of surface catalysis had been obtained for the 
TPS repair materials or damaged TPS material, ex- 
isting flight data from early Space Transportation 
System (STS) missions, which had highly catalytic 
C742 coatings on selected HRSI tiles, were used to 
obtain catalytic heating factors that could be ap- 
plied to any region on the Orbiter. 11 These C742 
coatings were placed on two selected tiles along the 
Orbiter centerline for STS-2 and STS-3. STS-5 con- 
tained C742 coated tiles at six centerline locations 
and two wing locations as indicated in Figure l. 2, 3,4 

During these early Shuttle flights, a unique phe- 
nomenon occurred when the flow passed from the 
partially catalytic RCG Coated tiles to tiles with the Figure 1. C742 catalytic tile locations on 
C742 coating. The heating rates on the highly cat- STS-2, -3, and -5. 
alytic tiles increased sharply above the partially cat- 
alytic tile values and even increased above the fully catalytic or equilibrium flow predicted values. The C742 
coated-tile heating rates then drop gradually across these tiles and drop sharply again as the flow then passes 
over the next partially catalytic tile. Predicted values using the BLIMPK code for Flights STS-2, -3, and 
-5 have been documented that show this overshoot phenomenon. 10: 12, 13, 14, 15 The overshoot is caused by a 
sudden drop in oxygen atom (O) concentration at the edge of the highly catalytic surface, with a consequent 
sudden increase in the oxygen molecule (O 2 ) concentration. This causes a sudden release of dissociation 
energy and increased heating to the edge of the highly catalytic tile. 

This phenomenon is not limited to tiles, but is documented for flow over the Orbiter RCC nose cap in 
which the SEADS pressure ports were installed. These ports were made of highly catalytic columbium metal. 
The sudden increase of chemical heating rate from the partially catalytic RCC surface as computed by the 
BLIMPK code over the columbium ports was explained by the sudden drop in atomic oxygen at the wall. 
After the flow passed the port, the oxygen atom concentration increased back to its nominal level over the 
RCC, and the heating was reduced to the lower catalytic value. The heating rate then alternately increased 
and decreased over each port in a saw-tooth fashion. 9 

The objective of this first catalytic heating study was to use the STS OEX flight data to determine 
catalytic heating factors by comparing the heating rates of HRSI tiles with the highly catalytic C742 coating 
to heating rates of the partially catalytic tiles with the original RCG coating. The study was based on the 
Development Flight Instrumentation data, 2,3,4 which contained the measurements from STS-2, -3, and -5. 

Recommended catalytic heating factors were defined for the Orbiter nose cap and chin panel region and 
the wing leading edge based on the catalytic factors derived from this study of STS-2, -3 and -5 OEX flight 
data. 16 These factors were later used to perform thermal analyses on the Orbiter RCC panels by the Boeing 
thermal group that supports the Shuttle Orbiter LESS/RCC NASA System Engineer, Dr. D. M. Curry, of 
NASA JSC ES3. 17 
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B. Arc-Jet Data 


A search was made for previous (historical) data at the JSC arc-jet facility as well as for more recent 
data obtained in the past year. This section describes some of the data which was converted into catalytic 
factors for a range of enthalpies. 

An assessment of wall temperature was made between coated and uncoated RCC test articles tested 
downstream of the exit of a conical nozzle in the JSC arc-jet Test Position 2 (TP2). 18 This test data was 
generated in 1990-91 and was reported for uncoated RCC 19 and for coated RCC. 20,21 Unfortunately, the 
test conditions were not the same between the uncoated and coated test articles. However, the data for the 
coated articles can be extrapolated back to lower values of the heating parameter given below 18 
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(. H c i - H w ) 


(1) 


Using this data, the measured surface temperatures were obtained from a laser pyrometer at two condi- 
tions of stagnation enthalpy for uncoated RCC and coated RCC. Catalytic (or bump) factors of 1.50-1.55 
were calculated 22 from the following equation, assuming the uncoated RCC and coated RCC emissivity were 
equal: 
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Recent tests have been performed at the JSC arc- 
jet facility in TP2 for uncoated or damaged RCC at 
higher enthalpies. The Pyrovision Infrared camera 3, 
was used to obtain the distribution of surface tem- 
peratures on the test article as a function of time 
for these tests. An example of a temperature dis- 
tribution on an RCC test article taken with the 
infrared camera during an arc-jet test is shown in 
Figure 2. In this case a fully catalytic damage site 
with exposed carbon substrate is located at the cen- 
ter of the test article. The results of these tests 
have not been published yet; however, preliminary 
results have been obtained 11 . Because of the range 
of calculated stagnation enthalpies (from 9,500 to 
13,000 Btu/lb) and the variation of surface temper- 
ature during the tests, it was not possible to obtain 
a precise value of enthalpy for each catalytic factor. 

An extreme case of damaged (uncoated) RCC was 
performed recently at TP2 in which the calculated 

stagnation enthalpy was about 20,000 Btu/lb - significantly higher than any flight value for orbital return. 
The range of catalytic heating factors from evaluation of arc-jet data was determined to be 1.62-1. 99. 22 



Figure 2. Temperature distribution on RCC test 
article from infrared camera during arc-jet test. 


C. Catalytic Recombination Coefficients 

In order to analyze the catalytic heating effects analytically, the recombination coefficients for the differ- 
ent surface materials must be defined and modeled correctly. Chemical heating occurs during atmospheric 
re-entry when the temperature of the air becomes large enough to excite and dissociate the air molecules 
into their respective constituents, namely oxygen and nitrogen. The catalytic heating depends on the re- 
combination rates of the oxygen and nitrogen atoms of the dissociated air at the wall and their respective 
heats of formation. If the wall is fully catalytic, all of the atoms will recombine to become molecules at the 
wall, releasing the heat of dissociation and increasing the heating rate to the surface. For this case, the flow 
field in the vicinity of the surface will essentially be in chemical equilibrium. On the other hand, if the wall 
is non-catalytic to recombination, none of the atoms recombine to become molecules. In this case, the flow 

a M9104 Mikron Pyrovision infrared camera (imaging pyrometer) manufactured by Mikron Infrared, Inc. of Oakland, NJ 

''Information transmitted from A. Rodriguez and R. Lewis/NASA-JSC to W. Rochelle/ESCG-Jacobs, 2005 
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near the surface is essentially frozen at the particular atomic mass concentrations present at the edge of the 
boundary layer. This produces the lowest heating rate to the surface. 

In reality, the surface of the HRSI tiles is neither non-catalytic nor fully catalytic, but partially catalytic. 
Hence, some of the atoms recombine at the wall, with resulting heating rates ranging between those expe- 
rienced for a non-catalytic wall and a fully catalytic wall. The glassy outer RCG coating of the Orbiter 
HRSI tiles tends to have low catalycity compared to metallic surfaces such as copper or columbium which 
are nearly fully catalytic. 

The recombination rates at the wall can be derived from the normal mass flux term in the Navier-Stokes 
equations for a gas in chemical non-equilibrium . 23 Early analyses used a catalytic recombination speed, k w , 
on the order of 700 cm/sec for both oxygen and nitrogen recombination at the surface. However, later studies 
generated separate recombination rates for oxygen and nitrogen recombination based on experimental data. 
These rates can be determined as a function of wall temperature from data in arc-jet facilities for RCG 
coatings on HRSI tiles . 6,24 


Temperature (K) 




(a) Recombination coefficients vs. temperature. 


(b) Comparison between two and four-equation RCG 
recombination models. 


Figure 3. Recombination coefficients and impact on heat transfer. 


General expressions for recombination rates developed for RCG-coatecl tiles 6 have been used in the litera- 
ture for the predictions of surface heating rates for several years. These expressions for oxygen recombination 
rate coefficients, 70 , for two ranges of surface temperature are 

7o = 40e _1144o ^ T ” 1435K <T W < 1580K (3) 

= 39 x 10 _9 e 2141 °' /T ™ 1580K <T W < 1845K (4) 

The nitrogen recombination rate coefficients, 7 n, for two ranges of surface temperature are 

7n = 6.1 x 10- 2 e" 248O/T ” 1410K <T W < 1640K (5) 

= 6.1 x lCT 4 e 5090/T ” 1640K < T w < 1905K (6) 

Combined, these expressions make up the two-equation temperature-dependent model for the partially cat- 
alytic recombination coefficients for the HRSI RCG-coated tiles. The two-equation temperature-dependent 
nitrogen and oxygen recombination coefficients for RCG , 6 C742, 7 and RCC (C-CAT ), 25 as well as the two- 
equation temperature-dependent air recombination coefficient model for RCG 25 are plotted vs. temperature 
in Figure 3(a). Although a four-equation model for RCG partially catalytic recombination coefficients is avail- 
able, the two-equation model was used in this analysis for consistency with previous analyses. A comparison 
between the four-equation and two-equation models for RCG recombination coefficients was performed to 
understand the differences in the two models. An example plot for this comparison is shown in Figure 3(b) 
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for a streamline out to the Orbiter wing leading edge. The results of the comparison revealed that the 
difference between the models is insignificant. 

A large database of temperature-varying recombination rates was generated in arc-jets and side-arm re- 
actors at NASA Ames Research Center several years ago. 25 Equations (3-6) as well as other equations in 
this database have been used in boundary layer programs such as the BLIMPK code, 5 and in CFD pro- 
grams such as the Langley Aerothermodynamic Upwind Relaxation Algorithm 26 (LAURA) and the General 
Aerodynamic Simulation Program (GASP) code 27 to predict heating rates on partially catalytic surfaces. 

D. Standard Two-Layer Methodology 

The general two-layer methodology is an inviscid/ viscous analysis consisting of an inviscid solution for 
the free-stream flow to the edge of the boundary layer and a viscous boundary layer analytical procedure 
to calculate the flow properties and heating effects from the edge of the boundary layer to the surface. 
The boundary layer analysis for this catalytic heating study was conducted with the BLIMPK code. 5 In 
past analyses with the BLIMPK code, the streamline inputs needed from the inviscid solution were usually 
extracted from inviscid Euler solutions. In a preliminary analysis, 28 the BLIMPK code was used with the 
streamline inputs extracted from an existing inviscid Euler solution obtained from the Inviscid Equilibrium 
Code in 3-Dimensions (IEC3D). 29,30 The streamlines were extracted along the Orbiter centerline for Mach 12 
and 18 in the STS-2 trajectory with the Axisymmetric Analog Two-Layer Convective Heating (AA2LCH) 
procedure. 31 Figure 4 shows the flight data compared to the BLIMPK analysis for partially catalytic RCG 
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(a) Heating rate ratio for Orbiter centerline. (b) Magnified. 

Figure 4. Standard BLIMPK/IEC3D two- layer methodology compared to OEX flight data for Mach 18 in the 
STS-2 trajectory. 


with C742 tiles and for full equilibrium for Mach 18. The forward most catalytic location is magnified in 
Figure 4(b). Because the preliminary BLIMPK analysis verified the thermo-chemical model in the BLIMPK 
code could satisfactorily model the catalytic heating overshoot observed along the centerline in the OEX 
flight data, NASA JSC requested that further analyses be conducted to model the catalytic heating effects 
observed on the Orbiter in the STS-2, -3 and -5 trajectories. However, since inviscid Euler solutions have 
largely been supplanted by computational fluid dynamics (CFD) Navier-Stokes solutions, a new method for 
obtaining the streamline inputs to BLIMPK was needed. 

The catalytic heating analysis presented in this paper was performed with a modified two-layer approach 
rather than the general two-layer method. Similar to the general two-layer approach, the modified two-layer 
approach consists of the boundary layer solution with the BLIMPK code, but uses extracted streamline 
inputs from three-dimensional Navier-Stokes (viscous) solutions at the edge of the boundary layer instead 
of inviscid Euler solutions. The streamline inputs are extracted with a streamline mapping algorithm called 
SRFTRACE, 32 described in more detail in the next section. The SRFTRACE algorithm is similar in nature 
to the AA2LCH procedure used to extract streamline inputs from the inviscid Euler solutions and calculates 
the same streamline inputs needed for the BLIMPK code including streamline distance, metric coefficient, 
pressure distribution and entropy gradient. However, the SRFTRACE code is also capable of extracting 
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the boundary layer edge species mass fractions from chemical non-equilibrium Navier-Stokes solutions, a 
capability not present in the AA2LCH code. With these streamline inputs, a non-equilibrium, non-isentropic 
analysis with defined boundary layer edge species mass fractions can be performed with BLIMPK. 

E. Modified Two-Layer Methodology 

In order to perform the catalytic heating analysis with the modified two-layer approach, three analytical 
tools were utilized in the following order: the Data-Parallel Line Relaxation (DPLR) CFD code , 33 the 
SRFTRACE code, and the BLIMPK code. 

1. DPLR 

The DPLR CFD code 33 was used to compute laminar Navier-Stokes solution at discrete trajectory points 
(time, altitude, and Mach number) for the Orbiter. All of the streamlines analyzed for the BLIMPK catalytic 
heating analysis originate from DPLR solutions. DPLR employs a thermo-chemical model and uses the four- 
equation temperature-dependent model for RCG for the entire windward acreage. Orbiter grid models exist 
for a low resolution (1-inch spacing between grid points) and a high resolution (0.25-inch spacing between 
grid points) computational analysis . 34 It is possible to analyze localized catalytic heating analyses for the 
Orbiter with the DPLR code, but it is time and resource intensive. Because catalytic heating analyses 
may be needed at discrete points while the Orbiter is in-flight, the time needed to complete a localized 
catalytic heating analysis with DPLR can be restrictive. Additionally, this time requirement precludes the 
use of DPLR for rapid turnaround, large-scale parametric studies like those considered in the present work. 
However, DPLR laminar solutions for RCG partially catalytic tiles exist for a number of trajectory points 
from which streamline inputs can be extracted. All of the DPLR solutions are for a nominal re-entry angle 
of attack of 40° and range between Mach numbers of 6 and 25 for the STS-107 trajectory . 35 The simulations 
employed a 5-species chemical non-equilibrium, thermal equilibrium model. DPLR-calculated heating rates 
are compared with the BLIMPK-calculated heating rates in section III. All of the streamline inputs for 
BLIMPK were extracted from the DPLR Navier-Stokes solutions with the SRFTRACE code. 

2. SRFTRACE 

The SRFTRACE code was used to map 
the data needed from the three-dimensional 
Navier-Stokes simulations into the format re- 
quired for the BLIMPK code. The princi- 
pal data inputs required for BLIMPK are the 
boundary-layer edge properties, the spatial lo- 
cation of the surface points, and the streamline 
metric coefficient. These first two input sets 
are trivially obtained from the CFD volume 
solutions, but the computation of the stream- 
line metric coefficient warrants further expla- 
nation. 

Figure 5 shows a representative streamline 
on the lower surface of the Orbiter. The ref- 
erence coordinate system is denoted by the 
{x, y, z) axes, which is the coordinate system 
used in the CFD solver. The (£, r], Q stream- 
line coordinate system is introduced to provide „ „ „ 

' .Figure 5. Jtieterence and streamline coordinate systems, 

the mapping tor BLIMPK. In this curvilinear 

coordinate system the £-axis is tangential to 

the streamline and the £-axis is normal to the surface of the vehicle. The 77 -axis completes the system. 

An arbitrary location can be expressed in the reference Cartesian coordinates as 

r = xi + yj + zk (7) 
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Furthermore, points on the surface of the vehicle can also be defined as 


r = + 'Qt-n 


( 8 ) 


The unit vector e v may be expressed as 


dr 

drj 


dr 

dr] 


( 9 ) 


where the magnitude term on the left hand side of the equation denotes the metric coefficient required by 
BLIMPK. This quantity measures the effective radius of curvature of the vehicle surface in the direction 
normal to the streamline. 

Consider the identity 

~dx, . dy, , dz u (1Q) 
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The last term on the right-hand-side may be expanded using the following identity 


dQ = dxd(J 
dy dr] dx 


dydQ 

dy dy 


dzd 0 
drj dz 


Using (12) in (11) yields the following system of ordinary differential equations 


d f dr 


dt \dy J dy 


dr 


7:hr =^'VV 
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(13) 


which defines an initial- value problem for the metric coefficient distribution along a given streamline. At the 
stagnation point of the vehicle the nose radius is known and defines the initial value of the metric coefficient. 
An explicit time integration technique is used to compute the metric coefficient at downstream locations, 
similar to the approach outlined in Reference 31. 

This procedure is implemented in a C++ code called SRFTRACE. This tool uses the libMesh c open- 
source finite element library to read the surface and boundary-layer edge parameters, along with the surface 
grid. The surface grid is treated as a hybrid-element unstructured mesh, which removes any constraints 
on the topology of the original grid. The standard piecewise linear Lagrange basis functions are used for 
interpolation along the surface. The boundary-layer thickness and edge properties are located using the 
curvature method described in Reference 36. 

The user specifies a point of interest on the surface and a streamline is generated via explicit integration 
in the direction opposite of the boundary-layer edge velocity. This procedure terminates in the upstream 
direction when the stagnation point is located. The streamline may also be extended downstream some 
user-specified distance if desired. The resolution of the points along streamlines depends on the resolution 
of the Navier-Stokes solution. In this study, the highest-resolution grids currently available for Orbiter RTF 
reentry 34 simulations were used, providing a streamwise point distribution on the order of 0.25 inches. At 
each point along the streamline the tool interpolates a number of values from the CFD data including (i) 
surface pressure, (ii) species mass fractions at the edge of the boundary layer, (iii) total energy, (iv) surface 
normal, and (v) boundary layer edge velocity. The velocity gradient is also computed, as is entropy. In 
the current implementation, the entropy is computed assuming local equilibrium using the curve fits by 
Tannehill. 37 It would be preferable to interpolate the entropy directly from the CFD solutions as well, but 
this is not supported with the existing post-processing tools. The tool outputs a table including all the 
information needed by BLIMPK, and optionally can produce a Tecplot® d binary surface file including all 
the derived parameters for inspection. 

c http : / /libmesh . sourcef orge . net 

d http://www.tecplot .com 
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3. BLIMPK 


The bulk of the catalytic heating analysis, including calculation of the catalytic factors, was performed 
with the BLIMPK code, a FORTRAN based code developed by Aerotherm for NASA for use in aeroheating 
analyses during the Apollo program. The BLIMPK code was used in this analysis to calculate localized 
catalytic heating at points along centerline streamlines and streamlines out to the wing tile acreage and 
wing leading edge. The BLIMPK model, described in the latest User’s Guide, 5 incorporates temperature- 
dependent catalytic models and utilizes 15 nodes across the boundary layer. Subsequent modifications allow 
the catalytic model to be updated at each station along the streamline, as well as the ability to model 2,500 
streamline stations. Macros and scripts have been used to speed up the pre-processing of the BLIMPK input 
files for large parametric studies. 


Table 1. Freestream and stagnation properties for STS-5 and STS-107 Best-Estimate Trajectories. 


Flight 

Time 
From El 
(sec) 

Altitude 

(ft) 

Mach 

No. 

Stagnation 

Enthalpy 

(Btu/lbm) 

Stagnation 

Pressure 

(atm) 

STS-5 

810.0 

208,560 

18 

6,886 

0.054392 

STS-107 

375.8 

248,430 

25 

11,839 

0.019555 

STS-107 

602.4 

230,907 

23 

10,324 

0.034366 

STS-107 

916.3 

203,819 

18 

6,830 

0.050651 

STS-107 

993.1 

192,972 

16 

5,619 

0.082356 

STS-107 

1339.6 

130,740 

6 

1,016 

0.142274 


The BLIMPK code input was set up to model a chemically 
non-equilibrium boundary layer of 5 species with finite gas reac- 
tion rates for non-isentropic flow around a body. The inputs also 
included streamline distance, metric coefficient, pressure distri- 
bution, boundary layer edge species mass fractions, and entropy 
gradient. In addition, BLIMPK required stagnation enthalpy and 
pressure values which were extracted from the Best-Estimated 
Trajectory for STS-5 4 and STS-107, 35 the two trajectories used 
in this analysis. These stagnation values are provided in Table 1. 

F. Orbiter Locations Investigated 

One of the goals of this analysis was to determine how well 
non-equilibrium, non-isentropic BLIMPK results based on a mod- 
ified two-layer methodology would match the STS-2, -3 and - 
5 OEX flight data. Because the STS-5 mission had six near- 
centerline C742 catalytic tiles in addition to the two wing C742 
catalytic tiles, the STS-5 Best-Estimate Trajectory was utilized 
for a major part of the present analysis. In order to verify that the 
BLIMPK analysis was suitable for predicting the catalytic over- 
shoot phenomena encountered during Orbiter atmospheric reen- 
try, all of the centerline C742 catalytic sensors were investigated, 
as well as the two C742 coated wing tiles on STS-5. The thermo- 
couple locations for the centerline sensors and two wing sensors 
on the Orbiter used for this analysis are listed in Table 2. 2,3,4 
Validating the current approach with these flight data allowed the 
methodology to be applied to the Orbiter nose cap, lower fuselage 
acreage, wing acreage, and wing leading edge with confidence. 


Table 2. STS-5 data sensor locations. 


Sensor 

Label 

Xo 

(in) 

Yo 

(in) 

Zo 

(in) 

V09T9341A 

267.7 

0.0 

306.1 

V07T9452A 0 

361.5 

-5.2 

285.6 

V07T9463A 0 

450.4 

0.0 

278.6 

V07T9464A 0 

486.5 

-6.4 

277.4 

V09T9381A 

565.0 

-4.3 

275.6 

V07T9468A 0 

620.1 

0.0 

274.4 

V07T9471A 0 

755.9 

0.0 

271.4 

V09T9521A 

878.9 

-12.8 

268.8 

V07T9478A 0 

1002.0 

0.0 

266.0 

V07T9481A 

1129.3 

0.0 

263.2 

V07T9487A 

1265.0 

0.0 

260.2 

V07T9489A 

1392.4 

0.0 

262.8 

V07T9492A 

1511.3 

0.0 

274.8 

V07T9689A* 

1222.6 

280.3 

283.9 

V07T9692A* 

1294.1 

277.9 

282.0 


° Denotes C742 Coated Tile 
* Denotes C742 Coated Wing Tile 
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Figure 6. Surface pressure contours and streamlines to wing 
C742 tiles. RCC panel locations shown for reference. 


The STS-5 BLIMPK wing analysis was 
compared to the STS-5 C742 highly catalytic 
flight data measurements mainly to judge the 
ability of the modified two-layer methodology 
to predict catalytic heating jumps on stream- 
lines out to locations on the Orbiter wing. 
The two streamlines to the catalytic tile lo- 
cations on the underside of the left wing on 
STS-5 pass very close to the bow shock/wing 
shock interaction region on the wing leading 
edge (in the vicinity of RCC panels 8 and 9), 
as shown in Figure 6. In past analyses, there 
have been problems running BLIMPK with the 
adverse pressure gradient that occurs in the 
shock/shock interaction region. 38 This cur- 
rent analysis proves that BLIMPK can handle 
the large adverse pressure gradients which oc- 
cur along the streamline near the shock/shock 
interaction region on the wing leading edge, 
provided that the appropriate boundary-layer 
edge properties are used. One observation 
from the current analysis is that it is critical to account for entropy swallowing in the presence of adverse 
pressure gradients when using this modified two-layer methodology. 

The other purpose of this document is to present the 
predictions of the BLIMPK surface catalytic heating anal- 
yses for fully catalytic repair/damage locations at specific 
RCC sites including the nose, chin panel, and wing lead- 
ing edge. The present BLIMPK surface catalytic heating 
analysis was performed for two RCC nose cap body point 
locations - one on the windward side and one on the lee- 
ward side - and one chin panel location. Streamlines for 
three panels with six zones per panel for the STS-107 tra- 
jectory were investigated, each passing through a respec- 
tive body point on the wing leading edge of the Orbiter. 

A schematic of the wing leading edge panel showing the 
zone locations is shown in Figure 7(a). e 

Table 3 contains the body point locations on the nose, 
chin panel and wing leading edge for the RCC catalytic 
heating analysis. The body points for panels 7 and 9 in 
Table 3 are taken along a normal cut to wing leading edge. 

Figure 7(b) f shows the panel 9 normal cuts in detail with 
the body point locations labeled. The body points for 
panel 11 are taken along a planar cut in X-direction at 
constant value of Y at 60 percent semi-span. 

The RTF BLIMPK analysis for these cases was per- 
formed to determine fully catalytic heating bump factors 
at different locations on the Orbiter RCC nose, chin panel 
and wing leading edge for specific times (Mach numbers) 
in the STS-107 trajectory. The BLIMPK-predicted cat- 
alytic factors for the RCC critical regions were calcu- , T . , , , , , . , , 

° Note: n/A denotes a nose cap body point location 

lated 38 for comparison to the original catalytic factors 

derived from the OEX flight data sensors for the STS-2, -3 and -5 C742 highly catalytic tiles. 11 


Table 3. Body point locations. 


Panel 

Zone 

Body 

Point 

Xo 

(in) 

Yo 

(in) 

Zo 

(in) 

n/A 

n/A 

114 

255.3 

0.0 

312.3 

n/A 

n/A 

1000 

236.0 

0.0 

338.5 

n/A 

n/A 

1020 

267.1 

0.0 

310.0 

7 

1 

5115 

1067.5 

199.6 

280.7 

7 

2 

5113 

1059.3 

209.9 

284.7 

7 

3 

5105 

1052.4 

218.5 

289.8 

7 

4 

5100 

1051.2 

220.1 

293.8 

7 

5A 

5155 

1052.4 

218.5 

300.6 

7 

6 

5163 

1059.3 

209.9 

312.8 

9 

1 

5514 

1105.6 

246.5 

285.2 

9 

2 

5512 

1101.4 

250.6 

287.0 

9 

3 

5505 

1095.0 

256.5 

291.9 

9 

4 

5500 

1093.6 

257.6 

294.5 

9 

5A 

5553 

1094.3 

256.4 

300.6 

9 

6 

5563 

1102.8 

246.9 

311.5 

11 

1 

2609 

1149.6 

281.0 

289.7 

11 

2 

2607 

1143.9 

281.0 

290.5 

11 

3 

2602 

1126.0 

281.0 

290.5 

11 

4 

2600 

1117.0 

281.0 

296.6 

11 

5A 

4602 

1126.0 

281.0 

307.1 

11 

6 

4607 

1143.9 

281.0 

315.0 


e Information transmitted from M. Fields/Boeing Houston to W. Rochelle/ESCG- Jacobs, 2005 
information transmitted from M. Barnwell/Boeing Houston to W. Rochelle/ESCG- Jacobs, 2005 
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Panel 9 

55% Semi-span (Normal Cut! 




(a) Schematic of wing leading edge panel showing zone locations. 


(b) Schematic of panel 9 showing nor- 
mal cut on RCC wing leading edge and 
body point locations. 


Figure 7. RCC panel geometry and body point locations. 


G. Assumptions and Cases Analyzed 

Several defining assumptions were made in order to perform the BLIMPK catalytic surface heating 
analysis. These assumptions are listed below: 

1. Non-equilibrium, non-isentropic (Tannehill equilibrium entropy curve-fits 37 ) solution with species mass 
fractions defined at the boundary layer edge (5 species model: N, O, NO, N 2 , 0 2 ) 

2. NASA Ames Research Center two-equation temperature-dependent models for RCG, C742 and C-CAT 
(RCC) recombination coefficients from References 6, 7, and 25, respectively 

3. Temperature-dependent RCG emissivity used for entire analysis 7 

4. C742 catalytic recombination coefficients for each catalytic tile location (approx. 8 inches long) on STS- 
5 windward tile acreage analysis and fully catalytic recombination coefficients for each repair/damage 
location (approx. 6 inches long for nose and chin panel and 8 inches long for wing leading edge) for 
RCC analysis 

5. Streamline cases for STS-5 trajectory at Mach 18 for the C742 highly catalytic tile locations (listed in 
Table 2) 

6. Streamline cases for STS-107 trajectory at Mach numbers 25, 23, 18, 16, and 6 for body point locations 
in six zones at Panels 7, 9 and 11 (90 total cases) 

The STS-2, -3 and -5 OEX flight data sensor locations are shown in Table 2, the body point locations 
for the STS-107 analyses are listed in Table 3, and the relevant stagnation properties for each trajectory 
point investigated are shown in Table 1. Catalytic location length and recombination coefficient variance 
sensitivity studies were conducted for selected cases. 
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III. Results 


The following section presents the results of the BLIMPK catalytic heating analysis performed on stream- 
lines for selected points in the STS-5 (for comparison to the OEX flight data) and STS-107 trajectories. All 
of the heating rate results are presented in non-dimensional form as heating rate ratios based on a reference 
heating rate. The catalytic heating factors are derived in the same manner as they are in Reference 22, as 
the fully catalytic heating rate divided by the partially (nominal) catalytic heating rate. The centerline, 
nose cap and chin panel streamline analyses are evaluated first, and then the wing acreage and wing leading 
edge streamline analyses are evaluated. 


A. Lower Fuselage Centerline Region 

The catalytic heating effects from the six C742-coated tiles on the Orbiter windward centerline during the 
STS-5 mission were evaluated 28 to verify that the BLIMPK code could model the catalytic heating overshoot 
observed on the STS-5 OEX flight data. Figures 8(a) and 8(b) show a Mach 18 comparison of the STS-5 


O 

CO 

ir 


x 




Orbiter x/L Orbiter x/L 


(a) Heating rate ratio for Orbiter centerline. (b) Magnification. 

Figure 8. DPLR and BLIMPK centerline analysis compared to STS-5 flight data for Mach 18. 


OEX flight data and the BLIMPK calculated heating rate ratios for the centerline partially catalytic RCG 
tiles with highly catalytic C742-coating at sensor locations indicated in Table 2. 

The entire streamline calculation and all the available centerline flight data is presented in Figure 8(a), 
and a magnification of the upstream C742 tile locations is shown in Figure 8(b). The flight data for the 
C742 catalytic locations is for the center of the 8-inch tile, where the BLIMPK catalytic heating has dropped 
significantly from the peak value at the upstream edge of the C742-coated tile. It is at this center location 
of each BLIMPK catalytic jump that the comparison with the STS-5 OEX flight data was made. 

Figures 8(a) and 8(b) show that the BLIMPK calculated heating rates compare very well with the STS-5 
OEX flight data for the RCG and C742 coated tiles. The catalytic heating calculated by BLIMPK at the 
first C742 tile location slightly under predicts the STS-5 data value at that location (although it is well 
within the uncertainty of the data for this location at this point in the trajectory, 6.3% 4 ), but the rest of the 
downstream BLIMPK catalytic heating jumps overshoot the STS-5 flight data. The BLIMPK overshoot for 
the downstream catalytic locations is expected because BLIMPK does not model the transient properties 
of the chemical heating occurring during reentry. The BLIMPK thermochemical model does not consider 
recombination of the atoms upstream. As the dissociated atoms travel downstream, they are partially 
exhausted by each highly catalytic tile location. Therefore, there are fewer atoms to recombine downstream 
at the other highly catalytic locations. Because there is more than one localized highly catalytic location 
along the centerline, this effect can be observed in the STS-5 flight data. 
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B. Nose Cap/Chin Panel Region 

A BLIMPK catalytic heating analysis was 
performed along centerline streamlines for the 
STS-107 trajectory to predict catalytic heating 
effects on the critical nose cap and chin panel 
RCC region. 22 The analysis was conducted for 
leeward and windward streamlines for the nose 
RCC region of the Orbiter. The analysis was 
done with both C-CAT (RCC) and RCG par- 
tially catalytic models with fully catalytic re- 
pair/damage sites at body point 114 on the 
windward nose, body point 1000 on the lee- 
ward nose, and body point 1020 on the chin 
panel (cf. Table 3). 

Figure 9 shows the Mach 18 BLIMPK heat- 
ing rate ratio results of the windward stream- 
line analysis with the DPLR predicted heat- 
ing rate ratios included for comparison. The 
DPLR and BLIMPK calculated heating rate 
ratios compare very well for both the C-CAT 
and RCG analyses, but at the stagnation 
point, the BLIMPK C-CAT (RCC) result com- 
pares better with the DPLR result than the 
BLIMPK RCG result. The fully catalytic jumps from the C-CAT analysis are larger than the jumps from 
the RCG analysis. Because the partially catalytic heating rates from C-CAT are also lower than the RCG 
partially catalytic heating rates, the catalytic heating factors based on the C-CAT analysis are larger than 
the catalytic factors based on the RCG analysis. The BLIMPK C-CAT partially catalytic heating rate 
increases abruptly after the second fully catalytic jump because of the transition to the RCG HRSI tiles. 

The two fully catalytic jumps shown in Fig- 
ure 9 were examined more closely in the Mach 
25 windward centerline analysis. The convec- 
tive heating rate (equivalent to the heating 
rate ratios presented in this paper) transferred 
to the surface, q CO nv is the sum of a conduction 
term, q con d and a chemical term, q c hem : 

qconv — qcond A qchem ( 1 4 ) 

The convective heating rate contributions 
from the conductive and chemical heating rate 
terms are shown in Figure 10 for the BLIMPK 
Mach 25 windward centerline analysis. The 
conductive heating rate is the product of the 
thermal conductivity of the material and the 
temperature gradient at the surface. The 
chemical heating rate is computed from a com- 
bination of unequal thermal and multicompo- 
nent diffusion coefficients and species concen- 
tration gradients. 23 If thermal diffusion is ne- 
glected, Equation (14) can be rewritten 39 

qconv = ~ A ) + E ( 15 ) 



Figure 10. BLIMPK Mach 25 windward centerlines 
heating rate components. 



Figure 9. DPLR and BLIMPK windward centerline 
analysis for Mach 18 in the STS-107 trajectory. 
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where the normal mass flux of the i th species, M* , may be expressed as 


Mi = p^Djj 
3 



( 16 ) 



Figure 11. BLIMPK Mach 25 windward centerline species 
mass fractions. 


Recommended April 2004 
C-CAT - Leeward Nose 
C-CAT - Windward Nose 
C-CAT - Windward Chin Panel 
Recommended March 2005 
RCG - Leeward Nose 
RCG - Windward Nose 
RCG - Windward Chin Panel 
Arc Jet Test Data 


where the summation is over all species. 

The chemical (or catalytic) heating occurs 
during hypersonic reentry and causes increased 
heating at the surface of the Orbiter. The air 
molecules (mostly molecular nitrogen and oxy- 
gen, or N 2 and O 2 ) are highly excited by the 
hypersonic shockwave and dissociate into their 
atomic constituents, N and O, which may then 
react to form nitric oxide, NO. At the surface, 
the atoms may recombine to form molecules, 
releasing their dissociation energy in the form 
of heat which is then transferred into the sur- 
face. The rate of recombination determines 
the amount of chemical, or catalytic, heat- 
ing. Some surfaces are very prone to mas- 
sive recombination (highly catalytic materi- 
als), while some are deterrents to recombina- 
tion (partially catalytic materials). By track- 
ing the species mass fraction concentrations, 
the level of recombination, or catalycity can 
be observed. 

Figure 11 exhibits the five species mass 
fraction concentrations on the Orbiter wind- 
ward centerline at Mach 25 in the STS-107 
trajectory. Between 0.0 and 0.5 feet and 1.0 
and 1.7 feet, C-CAT recombination coefficients 
are being used for the surface heating calcula- 
tions. In these areas, the species mass fraction 
concentrations appear nominal for hypersonic 
reentry. Oxygen is almost fully dissociated, 
shown by the high amount of oxygen atoms 
compared to the number of oxygen molecules. 

There are more nitrogen molecules than ni- 
trogen atoms in the partially catalytic region 
because nitrogen has a higher excitation tem- 
perature, and therefore, does not dissociate as 
quickly as oxygen. At the fully catalytic dam- 
age sites, between 0.5 and 1.0 feet and 1.7 and 
2.2 feet, there are significant increases in the 
N 2 , O 2 and NO molecules, while there are also 
significant decreases in the amount of atomic 
nitrogen and oxygen. The increase in nitrogen 
and oxygen molecules shows that much more 
recombination is occurring, hence the higher 

heating rates that can be seen in Figure 10. After 2.2 feet, the species mass fraction concentrations appear 
nominal once again, but are now higher than they were for the other partially catalytic regions because the 
RCG recombination coefficients were used for the surface heating calculations at this point. 

The BLIMPK catalytic heating analysis along centerline streamlines for the STS-107 trajectory was 
performed to obtain catalytic heating factors on the critical nose cap and chin panel RCC region to be 
used in real-time thermal analyses for future Shuttle flights. 22 Catalytic heating factors were calculated for 
fully catalytic repair/damage on the nose cap and chin panel region with RCG and C-CAT (RCC) partially 
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Figure 12. Catalytic factors derived from Orbiter centerline 
RCC nose cap analysis. 
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catalytic heating rates at the trajectory points listed in Table 1 for STS-107. The RCG and C-CAT derived 
catalytic heating factors obtained from this analysis are plotted in Figure 12 along with the Recommended 
Catalytic Model for the Orbiter Nose Region proposed in April 2004 (based on the OEX flight data analysis 
in Reference 11) and in March 20 05 40 as an upper limit (based on the STS-107 analysis with BLIMPK). Arc- 
jet test data points (described earlier) are also shown for comparison. As mentioned before, the catalytic 
factors predicted by BLIMPK based on the C-CAT partially catalytic heating rates are higher than the 
factors based on the RCG partially catalytic heating rates. In general, the catalytic factors obtained on 
the leeward side of the nose cap (where partially catalytic heating rates are lower) tend to be more severe 
than the catalytic factors obtained for the windward side. The catalytic factors also tend to decrease with 
decreasing enthalpy and as the damage location moves further away from the stagnation point. 

C. Wing Acreage Region 

A major focus of this BLIMPK catalytic surface heating analysis was the Orbiter wing tile acreage for 
STS-5. One of the major objectives of this analysis was to prove BLIMPK could successfully handle the 
large adverse pressure gradient induced by the the Orbiter bow-shock/ wing leading edge-shock interaction 
and determine how well BLIMPK could model the catalytic heating overshoot observed in the flight data on 
the two C742 highly catalytic wing acreage tiles. 41 Streamlines passing through these sensors are shown in 
Figure 6. 
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(a) Heating rate ratio for entire streamline. (b) Magnification. 

Figure 13. DPLR and BLIMPK predicted heating rates for OEX sensor location V07T9689A at Mach 18 in 
the STS-5 trajectory. 


The convective heating rate ratio results for the Mach 18 streamline out to the V07T9689A sensor on 
the Orbiter during STS-5 are shown in Figure 13(a) with the C742 catalytic jump at a streamline distance 
of approximately 85 ft. Both the DPLR solution heating rates and the BLIMPK solution heating rates are 
shown for comparison. STS-2 and -3 had very similar RCG partially catalytic heating rates at the sensor 
location, but only the heating rate at the sensor location for STS-3 is shown. The DPLR and BLIMPK 
generated heating rates are very close, and the BLIMPK solution was able produce the same peak heating 
rate in the adverse pressure gradient near the SSI region. Figure 13(b) shows the magnification of the 
V07T9689A sensor location along the STS-5 Mach 18 streamline. The DPLR heating rate ratio at the 
sensor location is slightly lower than the STS-3 partially catalytic heating rate ratio, while the BLIMPK 
partially catalytic heating rate ratio is even lower. The BLIMPK predicted value at the center of the C742 
tile (location of the flight data measurement) comes very close to the STS-5 flight data value. The BLIMPK 
prediction could possibly be higher if more data points were included across the tile, but it is still within the 
uncertainty of the flight data measurement (6.3 percent 4 ). 

The convective heating rate ratio results for the Mach 18 streamline out to the V07T9692A sensor on 
the Orbiter during STS-5 are shown in Figure 14(a) with the C742 catalytic jump at a streamline distance 
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of approximately 90 ft. Both the DPLR solution heating rates and the BLIMPK solution heating rates are 
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(a) Heating rate ratio for entire streamline. (b) Magnification. 

Figure 14. DPLR and BLIMPK predicted heating rates for OEX sensor location V07T9692A at Mach 18 in 
the STS-5 trajectory. 


shown for comparison. Again, although STS-2 and -3 had very similar RCG partially catalytic heating rates 
at the sensor location, only the heating rate at the sensor location for STS-3 is shown. The DPLR and 
BLIMPK generated heating rates are very close, but there is a larger difference between the peak heating 
rates in the two solutions at the adverse pressure gradient near the SSI region. Figure 14(b) shows the 
magnification of the V07T9692A sensor location. The DPLR heating rate at the sensor location is in close 
agreement with the STS-3 partially catalytic heating rate, while the BLIMPK partially catalytic heating 
rate is even closer. The BLIMPK predicted value at the center of the C742 tile is greater than the STS-5 
flight data value even with the uncertainty of the flight data measurement (6.8 percent 4 ). 

Although the BLIMPK heating rate values for the V07T9689A sensor streamline are slightly lower than 
the flight data measurements at the sensor location, the BLIMPK derived highly catalytic to partially 
catalytic heating rate bump factors are higher. The BLIMPK calculated bump factor (1.39) is slightly 
higher than the factor computed from the STS-5 highly catalytic heating rate and the STS-3 partially 
catalytic heating rate (1.33) for the V07T9689A sensor location. The BLIMPK highly catalytic to partially 
catalytic heating rate bump factor for the V07T9692A sensor location is also slightly higher (1.47) than 
the value computed from the STS-5 highly catalytic heating rate and the STS-3 partially catalytic heating 
rate (1.4). 

D. Wing Leading Edge Region 

The purpose of the BLIMPK analysis 41 along the off-centerline (wing) streamlines was to obtain catalytic 
heating factors for certain STS-107 trajectory points at locations on the RCC wing leading edge which could 
be applied in a thermal math model. The catalytic heating factors obtained from the new analysis would 
supplement the factors for the wing presented in Reference 11. Although 90 streamline cases were analyzed, 
the BLIMPK code was able to achieve converged solutions for only 85 of those cases. All streamlines for 
Panel 7, Zone 6 were unable to reach convergence, possibly due to the very deep expansion around the wing 
leading edge that occurs at that location. 

Of the 85 converged solutions, the worst case heating rates and temperatures occur during Mach 23 in 
the STS-107 trajectory in the shock/shock interaction region on the wing leading edge at Panel 9, Zone 
3. Figure 15(a) shows the worst case heating to the RCC wing leading edge at Panel 9, Zone 3 for Mach 
23 in the STS-107 trajectory. The 8-inch fully catalytic damage location is shown at a streamline distance 
of approximately 75 feet. The BLIMPK predicted heating rate ratios are very close to those predicted by 
the DPLR CFD solution up until a distance of about 20 feet at which point the two solutions diverge to a 
maximum difference at the peak heating location on the RCC wing leading edge. The maximum difference in 
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the heating between the BLIMPK and DPLR solutions occurs at the higher Mach numbers (the worst-case is 
Mach 25). The difference drops as the Mach number (or enthalpy) decreases. This difference may be because 
the current analysis assumes local equilibrium when computing the entropy gradient along the streamline as 
discussed previously. This approach should be most accurate near the stagnation point of the vehicle where 
the boundary-layer edge Mach number is low, and would become increasingly invalid as the flow accelerates 
rapidly around the vehicle, violating the assumption of local equilibrium. 



(a) Heating rate ratio for entire wing streamline. (b) Magnification. 

Figure 15. DPLR and BLIMPK results for fully catalytic region at Orbiter RCC Panel 9, Zone 3 at Mach 23 
in the STS-107 trajectory. 


Figure 15(b) shows a magnification of the fully catalytic damage site where the asymptotic peak and 
the exponential decay of the catalytic overshoot can be seen over the entire catalytic location. The peak 
catalytic factors and the catalytic factors at the mid-point were obtained for each converged solution and 
are presented in Reference 41. The real in-flight catalytic heating factor at the leading edge of the damage 
location is in between these two values. Although the worst-case heating rates and temperatures occur on 


Panel 9, Zone 3 (windward), the worst-case catalytic 
and 6), as will be discussed below. 

Figure 16 displays the mid-point catalytic bump 
factors versus streamline distance in feet for all zones 
(for converged solutions) for Panel 11. Figure 16 
shows that Zones 3 and 4 are closer to the apex of 
the wing leading edge and have the lowest catalytic 
factors, while the more windward and leeward zones 
have higher catalytic factors with the highest bump 
factor at either Zone 5 A or 6. The peak and mid- 
point catalytic bump factors for all converged cases 
along with the location and trajectory information 
and the fully catalytic and partially catalytic heat- 
ing rates are contained in Reference 41. The Mach 
6 cases were not evaluated as extensively as all fully 
catalytic bump factors were calculated to be unity, 
consistent with the negligible dissociation levels for 
this low freestream enthalpy (approximately 1,000 
Btu/lbm, cf. Table 1). 


factors actually occur on the leeward zones (Zone 5A 



Figure 16. Mid-point catalytic factors for Panel 11. 


17 of 21 


American Institute of Aeronautics and Astronautics 


In Figure 17, the recommended catalytic 
heating factor model for the Orbiter wing re- 
gion presented in Reference 11 is compared 
to the peak and mid-point catalytic factors 
predicted by BLIMPK in the present analy- 
sis for Zone 3 at Panels 7, 9 and 11, respec- 
tively. The recommended model was derived 
from the STS-5 C742 highly catalytic measure- 
ments on left wing tile locations and the STS-2 
and -3 partially catalytic measurements at the 
same locations. The BLIMPK catalytic factor 
predictions for the two wing catalytic tiles on 
STS-5 match very well with the recommended 
model as do the peak catalytic factors at Pan- 
els 7 and 9 from the STS-107 trajectory. The 
mid-point catalytic factor predictions at all the 
panels in Zone 3 from the STS-107 trajectory 
are significantly lower as are the peak catalytic 
factors at Panel 11. A limited amount of arc- 
jet catalytic data® obtained on a model placed 
in a wedge holder (flow parallel to surface) is 

also presented in Figure 17 for comparison with the BLIMPK predictions and the flight correlations. This 
range of arc-jet data resulted in factors of 1.47-1.66 based on a surface temperature increase of between 300- 
400 °F, respectively, on the damaged RCC coating (assumed fully catalytic) compared to the undamaged 
RCC (partially catalytic). The analysis of these factors was the same as that reported in Reference, 22 by 
taking the ratio of the fourth power of the temperatures and assuming the surface emissivity was the same 
for damaged and undamaged RCC. An undamaged RCC surface temperature of 2500 °F was used, and a 
computed arc-jet centerline (stagnation) enthalpy of 13,300 Btu/lb was determined. 



Enthalpy (Btu/lbm) 

Figure 17. BLIMPK peak and mid-point wing leading edge 
catalytic factors compared to recommended model. 


E. Sensitivity Study 

As noted in Reference 22, the catalytic heating factor is sensitive to many different parameters, and it 
is important to understand how those parameters affect the chemical heating. For the STS-107 analysis, 
catalytic location lengths were varied to determine if the size of a catalytic location had any effect on the 
overall bump factor. The recombination coefficients were also changed from RCG to C-CAT (RCC with 
Type-A Sealant) to observe the effect on the catalytic heating factors in the STS-107 wing leading edge 
analysis, although in this analysis, the flow would traverse RCG coated tiles from the stagnation point, then 
the RCC wing leading edge, and then the catalytic damage site. In an actual flight case, a wing leading edge 
streamline would start at the RCC nose stagnation point, then flow over the RCG HRSI tiles, and then back 
over the RCC wing leading edge and over the catalytic damage. 

Figure 18(a) contains the heating rate ratio result for the sensitivity study on the length of the fully 
catalytic location for Panel 9, Zone 3 at Mach 18 in the STS-107 trajectory. Figure 18(a) reveals that 
changing the length of the fully catalytic damage site from 8 inches to 3 inches has no effect on the peak 
catalytic heating or bump factor. Because the length has changed, however, there is an effect on the mid- 
point catalytic heating rate and the corresponding bump factor. The new mid-point distance is only 1.5 
inches from the peak for the 3 inch catalytic region and has a bump factor of approximately 1.24, while the 
mid-point of the 8 inch catalytic region is 4 inches from the peak and has a bump factor of 1.14. 

Figure 18(b) shows the heating rate ratio result for the sensitivity study performed on varying the 
recombination coefficients between RCG and RCC (C-CAT) for Panel 9, Zone 3 at Mach 18 in the STS-107 
trajectory. Figure 18(b) reveals that the RCC recombination coefficients yield a lower partially catalytic 
heating rate while generating a higher fully catalytic heating rate at the damage site as seen in the STS-107 
RCC nose analysis in Figure 12. The higher fully catalytic heating rate and lower partially catalytic heating 
rate yield a larger catalytic bump factor (approximately 1.9 at the peak; 1.3 at the mid-point) than the RCG 
recombination coefficients (approximately 1.45 at the peak; 1.14 at the mid-point). 

g Information transmitted from A. Rodriguez and R. Lewis/NASA-JSC to W. Rochelle/ESCG- Jacobs, 2005 
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(a) Catalytic length sensitivity study results for 
STS-107 Mach 18 Panel 9, Zone 3. 


(b) Recombination coefficient sensitivity study results for 
STS-107 Mach 18 Panel 9, Zone 3. 


Figure 18. Results from sensitivity studies. 


IV. Conclusions 

This paper presented the methodology and results of a two-layer surface catalytic heating analysis on the 
Orbiter RCG windward tile acreage and critical RCC nose cap, chin panel and wing leading edge panels. 
Comparisons of the OEX STS-2, -3, and -5 flight data (with the highly catalytic C742 coatings) along the 
fuselage centerline and wing tile acreage were made with the BLIMPK heating predictions. The BLIMPK 
predicted catalytic heating jumps compared very closely to the catalytic heating jumps derived from the 
OEX flight data. This method was extended to the RCC regions to calculate heating jumps for localized 
fully catalytic locations at assumed points of repair or damage on the critical RCC. 

The highest catalytic heating factors calculated for the RCC regions were slightly over 2.0 for the nose 
cap leeward side and on the order of 1.5 for the windward side. Catalytic heating factors computed at RCC 
Panels 7, 9, and 11 at six zones per panel for five free-stream Mach nos. were lowest (on the order of 1.2) in 
regions where the heating rate was the highest (Zones 3, 4). The catalytic factors were the highest on the 
leeward side (Zones 5A, 6) and at Zones 1 and 2 on the windward side where heating rates were lower. 

The factors presented in this paper were recommended for use with potential repaired and damaged 
areas of RCC for real-time analysis for future flights in which in-flight repair would have to be made prior 
to reentry or flown as-is. 
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